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PLANETARY SYSTEMS AND REAL PLANETARY NEBULAE 
FROM PLANET DESTRUCTION NEAR WHITE DWARFS 


Ealeal Bear 1 and Noam Soker 1 

ABSTRACT 

We suggest that tidal destruction of Earth-like and icy planets near a white 
dwarf (WD) might lead to the formation of one or more low-mass — Earth¬ 
like and lighter — planets in tight orbits around the WD. The formation of the 
new WD planetary system starts with a tidal break-up of the parent planet to 
planetesimals near the tidal radius of about li?. 0 . Internal stress forces keep 
the planetesimal from further tidal break-up when their radius is less than about 
100 km. We speculate that the planetesimals then bind together to form new sub- 
Earth daughter-planets at a few solar radii around the WD. More massive planets 
that contain hydrogen supply the WD with fresh nuclear fuel to reincarnate its 
stellar-giant phase. Some of the hydrogen will be inflated in a large envelope. The 
envelope blows a wind to form a nebula that is later (after the entire envelope is 
lost) ionized by the hot WD. We term this glowing ionized nebula that originated 
from a planet a real planetary nebula (RPN). This preliminary study of daughter- 
planets from a planet (DPP) and the RPN scenarios are of speculative nature. 
More detail studies must follow to establish whether the suggested scenarios can 
indeed take place. 


1. INTRODUCTION 


The first white dwarf (WD) where IR excess has been observed G29—38 ( Zucke r man & Becklin 


1987) was spectr oscopically observed also to contain unexpected metals in its atmosphere 


(e.g., 

Reach et al 

200 

9). Today more than 35 circumstel 

ar dusty disks around WDs are 

known (e.e., von 

RddcI et a 

1. 2007: Brinkworth 

et al. 2009; 

Farilii et al. 2010a 

b, 2 

012; Dufour et al. 

2010; 

Melis et al. 

2010, 

2011 

Girven et al. 2012 

Kilic et al. 

2012; Xu & Jura 1 

>012 

: Berefors et al. 

2014; 

Rocchetto et al. 

2015) 

. The composition of these WDs, as deduced from spectra, is 


consistent with accretion from asteroids or comets, but not with accretion from the ISM fe.g. 
Sion et al. 1990 : Zuckerman et al. 2003 : Jura 2003 : Jura et al. 2006; Kilic fc Redheld 2007 


1 Department of Physics, Technion - Israel Institute of Technology, Haifa 32000 Israel 















































































2 


Gansicke et ah 

2007, 

2008; 

Farihi et al. 

2010a) 

. IR excess from WDs (e.g. Becklin et al. 


2005 

Kilic et a 

2005 

, 2006; 

Jura et al. 2 

007alb, 

2009a; Mullallv et ah 

2007; Kilic & Redhcld 

200tI 

Jural 2008 

Farihi et ah 

2009,2010a 

b,2011 

.a; Farihi 2011b; Girven et ah 2012; Brinkworth et ah 

2012' 

is commonly attributed to dusty disks formed from tidal destruction of small bodies 


near the WDs (e. 




Jura et al. 2009b: Farihi et ah 


Zuckerman fc Becklin 


2010a 


1987; 


b, 2011 al; 


Graham et ah 


Debes et ah 


1990 


20121 


Brinkworth et al. 2009 


Veras et al.ll2014al ). 


The orbits of planets and the sub-planet bodies, here after planetary system objects 
(PSO), are mainly stable throughout the stellar evolution, after their formation and after 


dynamical settling flTsiganis et ahll2005 : Levison et ah 2011). However, planetary migration, 
close encounters and binary companions can act to destabilise systems at latter times. Fur¬ 
thermore, a stellar binar y co mpanion can form extremely eccentric planetary orbits through 


the Kozai mechanism ((Dav i es et al.l 2014). It is though that after the massive envelope ejec¬ 


tion by stars on the asymptotic giant branch (AGB), some orbits become unstable during the 


time i 

Debes & Sigurc 

sson 

2002; 

Nordhaus et ah 

2010 

Bonsor et al 

2011 

Veras et ah 

2011, 

2013; 

Vovatzis et ah 

2013; 

Frewen & Hansen 

2014 

; I 

Vlustill et al. 

2014; 

Veras & Gansicke 


2015). For example, the orbits of surviving planets increase due to the mass-l oss 


m a manner 


that might lead to instabilities of the orbits of PSO (e.g., IVeras et al.l l2014al: iMustill et ah 


20141) . These PSO might later be t i dally disrupted and collide with the WD remnant of 
the AGB star (e.g., Soker et ah 2010 : Veras et ah 2014a : Di Stefano et al. 2015 ). Due to the 
vari ety o f composition and orbital parameters other scenarios are possible as well. For exam¬ 
ple, Veras et ah (2014a) investigate the tidal disruption of highly eccentric asteroids as they 
approach WD and estab lishes that spherical asteroids break up and form highly collisionless 
rings. De bes et al. (120121) show that mass-loss from the star during the post-main-sequence 
stage can cause resonances between planetesimals and a giant planet. 

The relevant PSO can be classified into three groups as follows. 

(1) Small bodies of a progenitor planetary system like asteroids, comets and small moons. 
When tidally destructed these lead to the W D pollution as discussed above (e.g., 

20141 : iFarihi et al.ll2014j: IWilson et al.l 120141 and references therein). 


Bergfors et ah 


(2) Earth-like and icy planets can be destructed to planetesimals. Some planetesimals will 
be directly accreted to the WD. In section [2] we study the tidal destruction to planetesimals 
and propose that some surviving planetesimals might form sub-Earth planets. We term this 
process daughter-planets from a planet (DPP) since the new planetary system that is formed 
is composed of the tidally destructed planet. 

(3) Gas giant planets that contain sufficient hydrogen mass of Mh, p > 1O~ 3 M 0 , set a nuclear 
burning on the WD. In section [3] we argue that the end product might be an expanding 
nebula around a hot WD. This we term a real planetary nebula (RPN), as the nebular gas 
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comes from a planet. 

While previous studies investigated in detail PSOs in the first group, here we concentrate 
on PSOs in the second and third groups. We do not model these scenarios in detail, but 
rather provide a series of arguments to speculate on the final outcomes of the tidal disruption 
of bodies larger than asteroids, i.e., planets. In sections [2] and |H] we discuss two possible, 
somewhat speculative, outcomes from the tidally destructed planet. The initial evolution 
of unstable orbits is common to the two cases, but the final outcome is determined by the 
composition and size the initial planet. Icy planets, or Earth-like planets which are tidally 
destructed will lead to DPP and Gas giant planets will lead to RPN. We summarize these 
points in section El 


2. PLANETESIMAL FORMATION 


2.1. Tidal destruction 


Although tidal destruction of planets near WDs has been discussed before (see Sec. [Q) , 
we turn to evaluate the size of surviving planetesimals, something that was omitted from 
previous studies. We study the process of tidal destruction when an Earth-like planet or an 
ice giant reaches the tidal radius around a WD. Gas giant planets that contain hydrogen 
are discussed in section [3j The tidal destruction ceases when internal solid forces are about 
equal to tidal forces. Planetesimals are formed from the solid part of the planet. In ice giant 
planets the icy envelope is of lower density and will be tidally stripped before the solid core 
is destructed. 


We consider below the final destruction of the parent planet. In reality the planet can 
be destructed during a number of periastron passages. For example, Guillochon et al. (2011) 
found planets of > 0.25Mj to survive the first passage in their assumed orbit, although it 
already lose mass at that first passage. Since we study the latest evolution of the material 
of the destructed planet, our proposed scenarios work also for planet destruction in several 
passages. 


The tidal radius for a metal core is given by 


R t ~ CR* 


p*\ 

Pp) 


1/3 


~ 0.65 


C 

L26 


M* 


0.5 M 0 


1/3 


Pp 


5 g cm 


-3 


- 1/3 


Rq, 


( 1 ) 


where R* and M* are the radius and mass of the WD, respectively, p* is the average density 
of the WD, and pp is the density of the planet's metal core. C is a constant that d epends 
on the assumption used in deriving the expression (for more details see Harris 1996). Harris 
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(1996) states the range of values for C to be 1.26 — 2.9, depending on the properties of the 
destructed object. The tidal radius R t (which ranges from 0.65 — 1.5i? 0 ) meaning is that 
when the center of the planet reaches this distance to the center of the WD, tidal forces start 
breaking-up the planet. 

Equation (|T0 holds as long as there are no other forces beside gravity. However, when the 
shredded pieces of the core decrease, at a certain size, internal forces can hold them against 
gravity, similar to the way mountains exist on Earth. The maximum size is determined 
by the equality of the stress resulting from gravity to an internal property of the material 
call ed tens ile strength. It is usually expressed as ultimate tensile strength (UTS), a s (e.g., 
Bochnhar dtll2004 and reference within). UTS is the maximum stress a material can withstand 
(before breaking apart) while being pulled. 

In our setting, the maximum ‘weight’ the UTS can hold is W s = a s A, where A ps 
is the cross section area of the planetesimal. The ‘weight’ due to the differential gravity of 
the WD is W d ~ p pl g d Rp h where p v \ is the density of the planetesimal and the differential 
gravitational acceleration exerted by the WD at the tidal radius is 

2GM*-R p i 


9d 


Rl 


( 2 ) 


Taking W s > W d we derive the condition on the radius of a spherical planetesimal that is 
held by the UTS 


cr« 


flpi < —• (3) 

Ppl9d 

For a future use we define from this equation the maximum differential acceleration that the 
UTS can withstand 


O s 


9cr — 


Pp\Rpl 


( 4 ) 


We note that fractures can occur not only at the UTS limit due to tension but also at 
the limit of the compressive strength due to crushing. The compressive strength is usually 
higher than the tensile strength but for elastic composition can cause distortion and fractions 
in the material. The tensile strength for rocks found in San Marcos, for example , range from 
125 —580M Pa, depending on the composition flHuirong-Anita fc Thomas 1120041 ) . Meteorites 
which are considered as a rubble pile material have lower tensile strengths that range from 
2 — 62M Pa (c ompressive strength can range between 20 — 450 M Pa; for more details 


see 


Popova et ah 1201 1). Our core composition is unknown. 


We scale the UTS with cr s = 


100 M Pa and the planetesimal density to be similar to that of the Earth core ([Pepe et al. 
2013). 


Substituting equation 02} into equation (J3J) we derive the maximum planetesimal size 
that can be held by the UTS as a function of the tidal radius 
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R 


pi 


<T S R'l 

Ppi 2GM* 


0.5 


100 


M* 


0.5 M 0 


Ppi 


5 g s 


-i 


x 


CTs 


1 x 10 9 dyn cm _2 y \0.65i?(7 


Rf. 


km. (5) 


The simple calculation leading to equation (£>]) does not take into account the self gravity 
of the planetesimal. To do so we introduce to the equation the acceleration due to the gravity 
of the planetesimal: 

_ ^Gpp\R'i\ 

9pl ~ TRj ‘ 

The condition to hold the planetesimal against tidal break-up reads 


9d < 9p\ T Qcr- 


( 7 ) 


Substituting equations (J2J) , (3} and (E|) into equation (J7J), gives 


2 GAR 
R*t 


47T 

< y G Ppl + 


cr s 

Ppi-Rpi 


( 8 ) 


Setting a s = 0 yields equation CD Equation dHJ) sets an upper limit on the size of the solid 
bodies, planetesimal to Earth-like planets, to survive. In figure [Tj we plot this limit (solid 
blue line) on the planetesimal/planet size (radius) as a function of its distance from the WD 
center. 


As can be seen from figure |T] the planetesimals survivability is a function of their radius 
and their orbital separation from the WD. Dashed arrows depict the evolution route discussed 
here for icy or Earth-like planets. First, solid planet, or the solid core of an icy planet, is 
tidally destructed when it approaches the WD (black vertical arrow) and crosses the tidal 
radius at R t = 0.65 R@ in the setting used here. Planetesimals are formed around R t with 
decreasing size as tidal forces continue to act. Eventually, below the planetesimal radius 
given by equation (jHJ) and drawn in solid-blue line, UTS holds the planetesimal intact. Some 
planetesimals fall-in and accreted by the WD (red arrows). From the surviving planetesimals 
(green arrows), we propose here, a new sub-Earth planet might be formed. 


2.2. Evaporation 


Let us disc uss a spec ific system that might shed light on our proposed scenario and 
its limitations. Silvotti et al. ( 20141) claimed for the presence of three Earth-like planets 
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Fig. 1.— The orbital distance (r) vs. the planetesimal radius (R p i)- The figure represents the 
formation scenario of a DPP. The blue line depicts the tidal destruction radius of solid bodies 
with a density of p p i = 5 g cm -5 near a WD of mass M* = O.5M 0 as given by equation (J8]). 
Large bodies are held by their gravity, while small bodies by internal forces. In the region 
above the blue line the bodies are stable, while below the line they are destructed by the 
tidal force of the WD. The dashed lines present the evolution studied here. The vertical 
black-dashed line represents an Earth-like planet approaching the WD, until its distance is 
below the tidal destruction radius and it breaks-up to planetesimals. Some planetesimals 
are accreted (red lines), while some survive (green lines) and, we speculate here, merge to 
form lower mass sub-Earth planets. Some of the the surviving planetesimals are spread to 
larger distance of up to several solar radii. This is not drawn here. 
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orbiting the sdB star KIC 10001893. The respective orbital periods of the planets they have 
found are P\ = 5.273 hours, P 2 = 7.807 hours, and P 3 = 19.48 hours. This star is an extreme 
horizo ntal branch (EHB) star that burns helium in its core and contains a very little envelope 
mass ( Baran et al.1 2011; Ostensen et ah 2011). Spectroscopically, EHB stars are classified 
as hot subdwarf (sdB or sdO) stars. In order to become an sdB star, the red giant branch 
(RGB) stellar progenitor must lose most of its envelope. These stars retain a hydrogen 
envelope of less than about 0.001 M q flBaran et al. 2011 ). The small orbital separations of 
the planets from the sdB imply that the system went through a common envelope (CE) 
phase, or that the planets are second generation p lanets formed i n a p ost-CE phase with a 
star (e.g., |Peretsll 2010 l; iTutukov fe Fedorovall 2012 l) . ISilvotti et al.l (120141 ) assumes a standard 
sdB mass of M sdB = O.47M 0 , \ogg ~ 5.35, T sdB = 27500 K and i? sdB ~ O.24i? 0 . Therefore, 
the estimated luminosity of KIC 10001893 is L sc j B ~ 3OL 0 . 


Stellar binary companions (e.g., lHan et al.l 120021.120071) and planets flSokerl 119981 ) can 


perturb the envelope of RGB stars, mainly via a CE evolution, and lead to the ejection of 


most of the envelope (for a single star scenario see, e.g., Yi 20081) . The envelope ejection 
leads to the formation of an EHB star. If the claimed planets around KIC 10001893 are 
real, then there are two possible scenarios for the system formation. ( 1 ) A stellar companion 


planets formed. The se ar e second-generation planets (e.g. 


Beer et al. 

2004; 

Perets 

2010 


Tutukov & Fedorova 2012 ). (2) One or more planets entered the envelope of the RGB 


progenitor and removed its envelope to form the EHB star (Soker 


(Soker 

19981). Only the last 

e.g, 

Bear & Soker 

2011 b 

). It 


was an icy planet or li ghter. The Earth-like planets were formed from the tidal destruction 


of this last planet ( Bear etah 20 111) . This second possibility is according to the dynamical 


scenario proposed here. However, there is one limitation that prevents our proposed scenario 
to occur around very hot stars. 

Since these planets orbit an sdB star whose luminosity is higher than that of an old 
WD, their evaporation should be rapid. More significant, the planetesimals that are formed 
from the tidally destructed original planet are evaporated on a much shorter time scale than 
planets do. Below we present a simple calculation that shows that planetesimals will be 
evaporated in less than a year, shorter than the expected time for planetesimals to from a 
planet. We her by relay on two estima tions of evaporation time sc ales bv IStone et al.l (120151 ) 
for comets and by l Owen & W ul (120131 ) and Bear & Soker (I2011al) for planets: 


Stone et al. (2015) in their equation 11 estimate the timescale of a comet (that is com- 





































































posed of ice or rock) to completely sublimate as 


^ev 


167 r R c Qc,vP c r 2 


0.01 


J WD 


R 


pi 


Qc, 


V 


100 km J \ 3 x 10 10 erg g 


-l 


1 R& 


-^sdB 

30 L, 


yr, (9) 


where R c is the comet radius and p c is the comet density which in our case are equal to 
.Rpi and p p \ respectively, Qcy is the latent heat of transformation (similar for both ice and 
silicates), r is the orbital distance and Lwd is the luminosity of the WD, which for KIC 
10001893 is the luminosity of the sdB star. According to equation ([9]) based on Stone et ah 
(2015), the planetesimals will be completely sublimated within several days. 


We note that the Yarkovsky-O’Keefe-Radviesvki-Paddock (YORP) mechanism, that 
spins-up an asteroid to a break-up rotation velocity, wo rks o n_a much longer time scale for 
the planetesimals discussed here. Scaling the results of Veras et al. (l2014bl) to the sdB star 
discussed above for r = 1 R Q and R p \ = 100 km, we find the break up to occur at a time 
scale of > 100 yr. This is much longer than the evaporation time of these bodies. Indeed, 
the YORP effect is usually discussed with regards to small bodies with a size of 10 km or 
less ( Veras et al. 2014b). 


Owen fe Wu (120131) describe the evaporation rate of a planet caused by high energy 


photons (EUV and X-rays) near solar type stars as 

-hHE f T sdB 


M = 7 ] 


4 Gpr 2 


10 


3OL 0 


o 


-2 


Vr 

10~ 2 


g s 


-i 


( 10 ) 


where, The is the luminosity in EUV and X-ray. We defined r] R = 7/?/he, where 7 ] H e = 
Lne/-hsdB is the fraction of stellar radiation thal is_emitted in the EUV and X-ray, and 7] 
is the efficiency of t he pro cess discussed bv lOwen fc Wul (120131) . Scaling r/ based on figure 
12 of Owen & Wu ({2013 1. we have i] = 0.05 — 0.2. Since in our case we have an sdB 
star, we evaluate the number of energetic photons (A < 912A) from the entire spectrum as 
rj = 0.15. Therefore, r] R ~ 0.01. For this case a 100 km planetesimal will evaporate within 
less than 0.01 yr. We support t his app r oximation (Eq . JTUl) with the different evaporation 
rates presented in Bear & Soker ( 2011al ). IBear fc Sokerl ( 2011al ) calculate that at IRq from a 


0.5M 0 sdB star the sublimation rate exceeds 10 14 g s -1 (depending on the model considered). 
This lower limit of evaporation rate will completely evaporate planetesimals of 100 km in a 
timescale of a year. 


To sum marize. 


near a WD, 


Stone 


et al. 


(12015) took into account the evaporation rate of comets 


Owen & Wu (2013) calculated the evaporation rate of planets near solar type 


stars and Bear & Soker (|2011a) calculated the evaporation rate from a planet including 
recombination near an sdB star. All of these models suggest that the planetesimals formed 
at IRq near an sdB star will completely evaporate within a year. 
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The assembling process of a planet from planetesimals is complicated and the time scale 
depends on the size of the planetesimals, pr oper ties of the central star, the orbital distance 
and other parameters. Kenyon & Bromley ( 20141 ) calculated this timescale for different solar 
masses (0.1 M 0 , 0.3M 0 and 0.5 M & ) and evaluate the timescale for planetesimals of the order 
of 100 km to become an embryo of 1000 km to be T p i >em ~ 5 x 10 5 - 10 7 yr at 2.5 AU 
depending on the ini tial mass o f the disk and_the stellar mass. The Keplerian timescales for 
the cases studied by Kenyon & Bromley (2014) are 


Arbit ~ 6 


R 


2.5 AU 


M* 


0.5 M 0 


yr, 


( 11 ) 


We define the ratio 


c = — ,em 

^orbit 


io 5 - i(r 


( 12 ) 


The_ KejAerian timescale for the newly formed planetesimals in KIC 10001893 studied by 


Silvotti et al. ( 2014 ) is r or bit ~ 5 x 10 4 yr. Taking ( to be as in the equation (1T2|) we 


estimate the timescale of embryo formation to be T p i jem « 50 — 500 yr. According to our 
evaporation estimate the planetesimals will be completely evaporated within less than a year 
and hence no embryo can form. 

We conclude that in the case of KIC 10001893 evaporation prevents planetesimals to 
form sub-Earth pla nets. Furthermo re, the planets existence is debatable according to the 
evaporation rate is tone et al.lfioIHl . We note that a mag n etic field migh t shield Jupiter- 
like planets from extensive evaporation flBear fe Sokerll2012l) . ISilvotti et al.l (120141) suggested 
such magnetic fields inhibiting evaporation of the planets around KIC 10001893. It is not 
clear thought whether Earth-like planets will posses a strong enough magnetic field to slow¬ 
down evaporation. For planetesimals we do not expect that magnetic field will prevent 
evaporation. 


3. REAL PLANETARY NEBULA (RPN) FORMATION 


A real planetary nebula (RPN) is defined by us as an ionized nebula formed by mass 
lost from a giant star that was rejuvenated by a planet. Namely, the nebular gas was once 
part of a planet. We here propose, echoing ICorradi et al. (20151). that the collision or a tidal 
destruction of a gas giant planet with a WD can lead to a RPN. Plausible progenitor systems 
to the RPN are post-common envelope binaries (PCEBs) with plan etary syst ems around 
them. Properties of such sy stems are summarized bv IZorotovic fc Schreiberl (120131) and 
Schleicher fc Dreizler ( 2014 ). Schleicher fc Dreizler ( 20141) studied 12 PCEBs. elaborating 
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on NN Ser. A recent summary of the question on whether the planets are first-generation 
planets, i.e., were formed with the stellar binary system, or second generation planets (e.g, 


Vo 


scho w et al.l 1201411 . i .e.. fo rmed in the post - CE ph ase from a circumbinary disk, is given 


Bear fc Sokerl ( 2014 1 and Schleicher et al. ( 2015 1. We note that PCEBs are only one 


by 

possible formation route to RPNs 

ft is important to note that there is a question as to the real presence of planets in some 
of these systems. In particular, some sys tems, if real, seem to be dynamically unstable, e.g. 


HW vir, QS vir, and NSVS1425, (e.g., iHorner et al 


20121 


2013 


Wittenmver et al. 


2013 


Hinse et a,hi 120141 and in a recent paper Hardy et ah (2015) disproved the existence of a 
brown dwarf companion in V471 Tan. We here assume that such systems exist, and that 
instabilities can send some planets into the central binary system. The collision of planets 
with WDs should not be considered too speculative as we already mentioned in previous 
sections. 


Collision of a plan e t with a hot core of an AGB (or post-AGB) st ar was discussed before 


QHarpaz & Soker 


1994; 


Siess fc Liviol Il999l: iDe Marco fc Sokerl 120021) . There are differences 


between these studies and the RPN scenario. (1) In the above listed studies the contribution 
of the planet to the ejected mass was l ower t han half. The largest contribution of the 
planet mass is in the calculations of Harpaz & Soker (19941), where for four cases a planet 
of mass O.O1M 0 was destroyed inside an AGB envelope of mass O.O2M 0 . (2) The planet in 
these studied were destroyed near a hot core that still has a hydrogen burning in a shell. The 
injection of more hydrogen can change the burning properties, but the burning already exists. 
In the RPN scenario the hydrogen from the planet sets a new burning phase on the WD. 
Namely, in the RPN scenario the gas from the planet is fully responsible for the reformation 
of an AGB-like star. The relevant process from these studies is that the accretion on to the 
WD is via an accretion disk. This can lead to the launching of jets, and to extra mixing of 
the outer layers of the core. 


When a planet of mass M p ~ 1 — 10Mj, where Mj is the mass of Jupiter, approaches a 
WD, or an sdB star, it is tidally destroyed at radius R t as given by equation (P). First the 
hydrogen-rich envelope of the planet is removed, and then its metallic core. As discussed 
in Section [2] the metallic core is actually fragmented to many planetesimals. However, no 
planet will be formed as the hydrogen starts to burn on the WD. The high temperature 
evaporates the planetesimals. 


About half of the gas from the destructed planet might escape the system, but about 
half of the gas stays bound and forms an accretion disk. The accreted hydrogen gas starts 
to burn on the surface of the WD. The accretion time scale can be up to hundreds time the 
dynamical time scale at the tidally destruction radius, i.e., several weeks. This leads to an 
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accretion rate of M acc _ t > 0.01 Mq yr h A WD cannot accommodate such a rate, and a 


red giant envelope is inflated ( Nomoto et al.l 20071) . A post AGB star with a core mass of 
~ 0.6 Mq maintain a giant envelope (R > lOOi? 0 ) as long as the envelope mass is > 10” 3 M 0 
flSokerl 1 1 99211 . Aro un d a cold WD of ~ 0.6M^ even a lower mass of 10 4 M 0 is sufficient to 
inflate a giant envelo pe flHachisu et al. I1999T ). More massive WD require less mass to inflate 
an envelope flHachisu et ah 1999) . 


In cases where there is a PCEB that causes the planet-WD collision, there is now a 
binary system inside the rejuvenated AGB star. Namely, a CE structure. The nebular mass 
is much smaller than any of the two stars. This situation resembles in some respects nova 
ejection that covers its binary system progenitor. To learn about the outcome we turn to 
examine nebulae formed by novae. 


Many nova ejecta posses b ipolar (or elliptical) shapes fe.g.. lWoudt et al.ll2009l : iGhesneau et al 
2 0111. 2012; Shore et ah 2013). The expansion velocity of novae is typically high, v e > 

500 km s _1 as the hydrogen or helium burn in a thermonuclear run-away, and the mass is 
very low^ much below 0.001M Q . However, in some cases slower components are observed 
(e.g., Chesneau et al. 2011). In the systems we study here the accreted hydrogen is not 


degenerate and the mass is much larger. Hence, expansion velocities will be much lower. 

It is thought that the CE during the nova ejection sh apes the bipolar nebula with 


possible mass concentration in the equatorial plane (e.g., Livio et al. 1990; Parcscc et al. 


1995). It is quite possible, therefore, that the bipolar structures of RPNe might result from 
a CE structure formed during a planet destruction. The major differences between novae 
and RPNe are (i) the ejection speed in novae are much higher, and {ii) the ejecta mass in 
novae is much lower. But the basic formation of bipolar nebulae might be similar. 


The possibility of an RPN was already discussed by Corradi et ah (2015) as one possible 
explanation for the three PNs Abell 46, Abell 63, and Ou5. These have post-CE close binary 
systems, and very low nebular mass. The mass in each of these nebulae is much lower than 
the AGB envelope mass required to bring a wide binary system to a close orbit. Hence, 
Corradi et al.l (2015) raised the possibility that the CE phase of these binary systems took 
place a long time ago, and the present low-mass nebula is a remnant of a planet. The general 
structure of these three PNe is bipolar with an equatorial ring. 


In principle the RPN scenario does not require a close binary system. But it does 
require some instability in a planetary system to cause a planet collision with the WD. A 
Kozai-Lidov mechanism induced by a far-away binary companion is another possibility to 
cause a planet to collide with a WD. In the present study we do not consider the dynamical 
evolution that leads to planet-WD collision, and that is the same for the DPP and RPN 
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scenarios. 


4. SUMMARY 


We conducted a study of planets that are tidally destructed by WDs. Although rare, 
these types of transient events are expected to be found with future surveys, and be able to 
teach us about planetary systems around evolved stars. 

The destruction of asteroids was studied in detail in the past. In this paper we speculated 
on the consequences of tidal destruction of planets. We did not model these scenarios in 
detail. We rather suggested plausible scenarios to the final outcomes of tidal disruption of 
bodies larger than asteroids. Earth-like planets and cores of icy planets which are tidally 
destructed might lead to the formation of daughter-planets from a planet, the DPP scenario. 
The newly formed daughter-planets are low-mass planets (Earth-like planet or less). Gas 
giant planets that are tidally destructed might cause the formation of real planetary nebulae 
(RPNe). Although the outcomes are different and are the direct result of the composition 
and mass of the planet (i.e., the hydrogen amount) both outcomes result from the same 
dynamical evolution leading to a tidal destruction of a planet. 

When an Earth-like planet or an ice giant reaches the tidal radius around a WD, the 
metal core is tidally destructed and ~ 100 km planetesimals are formed. The planetesimals 
survivability is a function of their radius and their orbital separation from the WD. Some 
planetesimals will be accreted by the WD but some that survive might merge to form one 
to few new sub-Earth planets. For ice giant planets the icy envelope is of lower density and 
will be tidally stripped before the solid core is destructed. The solid core then evolves as an 
Earth-like planet and might lead to the formation of a planetary system. 


Gas giant planets contain hydrogen that complicates and enriches the outcome. Part of 
the hydrogen that is accreted onto the WD starts to burn. The now hot WD evaporates any 
surviving planetesimals from the tidal-destruction process. The tidal destruction process and 
the friction inside the accretion disk that is formed can supply hydrogen at a rate much higher 
than what the WD can accommodate. Part of the hydrogen then inflates a giant envelope. 
This rejuvenates giant star loses part of the envelope in a wind that forms a nebula. Later 
the central objects shrinks, heats up and ionizes the nebula, like in a planetary nebula. We 
term this nebula whose material comes from a planet, a real planetary nebula (RPN). This 
scenario was mentioned by Corradi et al. (2015) in their study of some very low-mass PNe. 


The collision or tidal destruction of a planet near stars beyond the main sequence 
can occur not only around cool WDs, but around hot sdB stars as well. In this case the 
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planetesimals will be evaporated due to the high luminosity of the hot sdB star and a planet 
is unlikely to form. 

We thank an anonymous referee for very detailed comments that substantially improved 
our manuscript. 
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